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Abstract. We devise a microscopic scattering approach to probe the excitation 
spectrum of a Bose-Einstein condensate. We show that the experimentally accessible 
scattering cross section exhibits universal Ericson fluctuations, with characteristic 
properties intimately related to the underlying classical field equations. 
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1. Introduction 



Bose-Einstein condensates (BEC) in optical lattices provide a versatile tool to address 
experimentally a variety of questions that emerge in diverse fields ranging from 
quantum information and many-body quantum phase transitions to solid-state transport 
and atomtronics. An important element of these studies is the development and 
implementation of methods which allow for an accurate measurement of the properties of 
the condensate. Among the most popular ones are time-of-flight and Bragg-spectroscopy 
[U [2j [3] techniques which result in the destruction of the BEC, whereas only few works 
consider a scattering setup that leaves the condensate intact Specifically, the main 
focus of the present literature is on photon-atom scattering, while only very recently 
matter-wave scattering was experimentally implemented to infer the spatial ordering of 
ultracold atomic crystals [5], and theoretically proposed to probe the Mott insulator 
to superfluid transition of the condensate ground state [6]. However, excited states are 
naturally populated in experiments which probe non-trivial BEC dynamics [7]. The 
rapidly emerging complexity of the many-body dynamics - manifest, e.g., in dynamical 
instabilities [3] - is a direct fingerprint of the complex underlying spectral structure, 
which is itself reflected in the - in general chaotic - classical limit of the Bose-Hubbard 
model (achieved in the limit of large particle numbers). It is therefore timely to explore 
possible experimental strategies to probe these spectral features, in a non-destructive 
manner. In our present contribution, we show how an inelastically scattered probe 
particle can reveal the state of a BEC target in the parameter regime of spectral chaos. 
Due to the inherent sensitivity of spectral cross sections under such conditions, a robust 
characterization requires a statistical approach, which can be further sharpened by 
semiclassical considerations. 



2. Model 

The scattering setup that we have in mind is shown in Fig. [TJ A probe particle with 
momentum k moves in a waveguide which is placed in the proximity of a BEC confined 
by an optical lattice. When the particle approaches the condensate, it interacts with the 
latter - much as the condensate particles between themselves - leading to an exchange 
of energy. The particle energy on exit from the waveguide defines the scattering cross 
section. The dynamics of the process is generated by 

Htot = H TB ® i + i <g> # BH + H int , (1) 

where 1 denotes the identity operator. In (fij, if B H represents the BEC target's Bose- 
Hubbard-Hamiltonian [8] 

U L 

H B}i = -J2 H^i ~ !) - K J2 [$kf * + h - a ] ( 2 ) 

i=l i 

of N interacting bosons on an L-site optical lattice, with the bosonic annihilation 
(creation) operators, and n« = b\bi the particle number at site i. U and K parameterize 
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the on-site interaction and the inter-site tunneling strength, respectively. In the 
macroscopic limit N — > oo (UN fixed), the dynamics of the condensate is well described 
by mean-field theory, i.e. the discrete Gross-Pitaevskii equation. In this limit, the 
quantum operators bf' are replaced by L complex amplitudes of a single-particle 
field. The Hamiltonian ^ then reads 

^E = ^ElA| 4 -^E[A*A + i + c.c.], (3) 

1=1 i 

where the A{ are time-dependent, and obey the canonical equations idAi/dt = 
dH GP /dA*. 

The waveguide in our scattering scheme is (for mathematical convenience) modeled 
by two semi-infinite tight-binding (TB) leads with hopping term J and lattice spacing 
a — 1. The particle's energy in the momentum eigenstate \k m ) of each lead is 
e m = — 2Jcos(fc m ), with corresponding velocity v m = 2Jsin(fc m ) [9]. These two leads 
are coupled with strength J to the central site j = 0, which is closest to the BEC. 
Jo thus controls the effective coupling of the projectile-target interaction region to the 
asymptotically free states of the lead. The projectile Hamiltonian is: 



TB 



jy-i,o i=-i,o 



h.c. , (4) 



with Cj the annihilation (creation) operators of the probe particle at site j of the TB 
lead. 

Finally, the probe-target interaction Hi nt is assumed to be of similar type (i.e. short 
range) as the bosonic inter-particle interaction in the condensate: 

H int = ol c c ® Wi . (5) 

For non-vanishing tunneling coupling K, H[ nt induces transitions between different 
eigenmodes of the condensate, what renders the scattering process inelastic. In this 
sense, a > controls the inelasticity. As the setup does not allow particle exchange, 
target bosons and probe particle are distinguishable at all times. Thus, the latter 
can be chosen at the convenience of the experimentalist to implement the interaction 
described by ([5|. In the macroscopic limit, the interaction Hamiltonian becomes time 
dependent, and is given by Hmt/N = a\Ai\ 2 c\cq. 



3. Scattering matrix 

Given the total Hamiltonian and the asymptotic freedom of the probe particle, 
we can now define the scattering matrix of our problem, as the fundamental building 
block for our subsequent observations: For the BEC initially prepared in an energy 
eigenstate \E m ), and the probe particle injected with an energy e m , the total system 
energy is £ = E m + e m . The open channels (modes) of the scattering process are then 
determined by energy conservation and characterized by the kinetic energy e n = £ — E n 
of the outgoing probe particle. The transmission block of the scattering matrix can be 
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Figure 1. (color online). Scattering setup: The probe particle is injected into 
a wave guide, and locally exchanges energy with a BEC confined by a three-site 
optical potential, in the contact region between wave guide and site one of the lattice. 
The inelastic scattering cross section measured on the wave guide's exit bears direct 
information on the state of the BEC. 



derived from the Green's function of a particle at site j = 0, with two semi-infinite leads 
attached, and reads [|J 

[St](£) = V-v- — * 7 — ——, Ve, (6) 

(1 - 7j|£ - H B h\ - emi + ijv 
where 7 = (Jo/ J) 2 , and v is the velocity operator. In the eigenbasis of the BH 
Hamiltonian, [HBH\nm = E n 8 nm and v nm = v n 5 nm are diagonal matrices, while 
Q nm = (E n \hi\E m ) is not. For 7 = 1, Eq. ^ coincides with the S"-matrix for inelastic 
electronic scattering derived in [10]. In our setup, 7 < 1 can be regarded as a potential 
barrier that reduces the coupling between the leads and the scattering region (i.e. for 
7 = the latter is isolated and the probe particle is perfectly reflected). As 7 is increased 
from zero, one observes a crossover from a regime of well-resolved, narrow resonances to 
a regime of overlapping resonances. Whereas the former regime would in principle allow 
us to directly infer the BH spectrum, it puts very hard demands on the experimental 
precision as far as preparation and measurement are concerned. Instead, we find that 
in the latter regime one obtains an experimentally robust, fluctuating scattering signal 
that bears information on the target's spectrum (around 7 = 0.001). This information 
is lost when 7 becomes very large compared to the mean resonance spacing such that 



I The derivation of this formula in the present work is a tight-binding generalization of the continuum 
limit that has been analyzed in [TU]. Details on the derivation can be found in 
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Figure 2. (color online). Transmission T m , averaged over the channels m = 400 — 430 
in the chaotic regime, versus the total energy £ in units of K, for 7 = 0.001, N = 38 
bosons, and a — 5. OX. 

these characteristic fluctuations are washed out. 
4. Chaotic scattering 

In our present contribution, we will investigate the properties of a probe particle 
scattering on a BEC that is described by the Hamiltonian In contrast to 

Ref. [6], we are not concerned with the well-understood ground-state properties and 
the associated superfluid to Mott-insulator transition, which takes place at rather large 
interaction strengths U. Instead, we focus on the complex properties of excited states 
at intermediate interaction strengths. 

Namely, for L > 2 and values (3 < u < 12) of the control parameter u = UN/2K, 
the classical Hamiltonian Hgp ftty generates chaotic dynamics [121 [THJ [HI HHJ fT6] . 
Quantum manifestations thereof were investigated in a series of publications with 
emphasis on the statistical properties of the energy spectra [TFJ [TEl [19]. For a two- 
site lattice (L = 2) [20J , inelastic scattering revealed immediate fingerprints of the fully 
integrable mean-field dynamics. Here, we consider a three-site BH Hamiltonian (L = 3), 
intermediate values of the control parameter around u — 5, what corresponds to the 
maximally chaotic regime [12], large filling factors of the lattice (i.e. N = 38 bosons), 
and a condensate that is initially prepared in an energy eigenstate \E m ) in the bulk of 
the BH spectrum, where the target dynamics is predominantly chaotic. The presence of 
chaos qualitatively alters the resulting physics: In stark contrast to the integrable case 
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Figure 3. (color online). Left: Logarithmically color-coded snapshot of the underlying 
interaction matrix Q n rn, for u — 5. Right: Distribution P(Q nrn ) of the unfolded off- 
diagonal interaction matrix elements Q nm taken from a sub-matrix in the bulk of 
Qnm (indicated by the white box in the left panel; the corresponding indices run from 
n — m — 360 to n — m — 500). The distribution P{Q nm ) nicely resembles the dashed 
red line which represents the Gaussian distribution with unit variance N(0, 1). 



L = 2, the scattering quantities will strongly fluctuate (depending, e.g., on the probe 
particle's energy), what requires a statistical analysis of the scattering signal. On the 
other hand, chaos is expected to yield universal behavior, i.e., the obtained results should 
not depend on the details of the target system. As for the other scattering parameters, 
we set J = E max — E min = 211K, equal to the (numerically inferred) spectral width of 
the BH Hamiltonian, such that all scattering channels are open, and fix 7 = 0.001. 

How do the chaotic spectral properties of the BEC manifest in a scattering 
experiment as sketched in Fig. HI? The experimentally most easily accessible observable is 
the transmission T m {£) = J2 n \[^T]nm\ 2 - It denotes the probability that a probe particle 
with incoming energy e m exits the scattering area in any of the outgoing channels e n . 
In Fig. [2j we show T m versus the total energy S, averaged over 30 incoming channels 
e m in the middle of the spectrum, i.e. in the chaotic regime, where approximately 180 
outgoing channels contribute to T m : 

Strong fluctuations dominate the transmission signal - a first indicator for the 
complexity of the target. In order to understand how this complexity enters in the 
scattering process, we turn to the scattering matrix ^ . In the {\E n ) } basis, in which the 
latter is evaluated, the interaction operator hi becomes the only non-diagonal quantity 
on the r.h.s. of and thus gives rise to inelastic scattering processes. A closer 
inspection of the corresponding matrix Q nm = (E n \hi\E m ) (see Fig. |3^i)) shows that 
for intermediate energies corresponding to the center of the matrix (i.e., in the chaotic 
energy regime) the matrix is banded and its off-diagonal elements look rather erratically 
distributed. We analyze these elements, taken from the box indicated in Fig. |3^i), 
and consider their distribution P(Q nm ). Here, Q nm = Qnm/ a/ (\Qnm\ 2 ) represent the 



Inelastic chaotic scattering on a Bose-Einstein condensate 7 




Figure 4. (color online). Left: Integrated total inelastic cross section J d£{p^(S)}, 
for the same parameters as in Fig. [2j averaged over the channels m = 400 — 430 in 
the chaotic regime, versus the inelasticity parameter a in units of K. The integration 
runs over the entire energy axis. Middle: A representative inelastic cross section a nm 
versus the total energy £ in units of J, for the same parameters and a = b.QK. 
Right: Histogram of the normalized inelastic cross section P(o nm ), for fifteen different 
channels UAmm ( m — 401 — 415) in the chaotic regime and identical parameter values. 
The data perfectly match the dashed straight line exponential fit. 



off-diagonal matrix elements which are rescaled by their local standard deviation, the 
latter being calculated in a small, moving energy window. This "unfolding" is necessary 
to remove system-specific features and reveal universal properties of the interaction- 
operator hi (see, e.g. [12]). The resulting distribution is shown in Fig|3}D) together 
with a Gaussian of unit variance. The good agreement with the latter indicates that 
(neglecting higher-order correlations) the matrix elements are essentially independent 
random numbers, in perfect agreement with the predictions for quantum systems that 
possess a chaotic classical counterpart [2H [221 [23] . We conclude that St and thus all 
scattering quantities inherit their complexity from Q, since the latter represents the key 
ingredient in the scattering matrix. 

To gain insight in the role of the parameter a that controls the inelasticity induced 
by Q, we next consider the total inelastic scattering cross section 

p[^) = 2]T|[M£)UI 2 , (7) 

which essentially resembles T m , except for the direct processes. For a given value of a, 
we integrate over the energy axis, to obtain robust results, unaffected by the sensitive 
energy dependence of p^(S). Fig. |4^i) shows that J d£(p™(£)) m takes its maximal value 
for intermediate values of the inelasticity parameter a, while it vanishes in the limit of 
small and large a. In the former case, the probe particle is directly transmitted, since 
(g) with a ^ becomes diagonal, while it is directly reflected in the latter case - as 
evident from ^ with a IK in the denominator. Consequently, only for intermediate 
a-values around a ~ 5K can we expect to infer information on the condensate from the 
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probe particle's exit energy. 
5. Ericson fluctuations 

Beyond total cross sections there is nontrivial dynamical information encoded in the 
partial inelastic cross sections a nm {£) = |[SV(£)]nm| 2 , which quantify the probability 
for a transition from a state E m to a state E n of the target (or, equivalently, from an 
energy e m to an energy e n of the projectile). In Fig. |2]d) we show cr nm (£), for the same 
parameter values as the transmission in Fig. [2} We observe much stronger fluctuations 
than for the total transmission, what is simply due to the fact that the latter implies 
an additional effective averaging over many scattering channels. As we will show now, 
this sensitive dependence on the energy is an unambiguous trait of (universal) Ericson 
fluctuations, that were hitherto only reported in the context of nuclear [24] and atomic 
physics [25J, as well as in microwave experiments [26], and are also connected to universal 
conductance fluctuations in mesoscopic physics [2"?] |2"8"]. 

The rapid fluctuations of the cross section are due to interference effects between 
overlapping resonances: The scattering amplitudes [5r] nm can be represented by a 
many-resonance Breit-Wigner formula, where each individual term in the sum is a 
random variable whose stochastic properties originate from the statistically independent 
Gaussian distributed elements of the interaction matrix. Then, due to the central 
limit theorem, one expects that real and imaginary part of the scattering matrix 
elements are Gaussian distributed random numbers with zero mean. In other words, 
the Gaussian distribution of the interaction matrix elements gives rise to Gaussian 
distributed real and imaginary part of [SV] nm . This results in an exponential distribution 
[21] P(a nm ) = exp[—a nm ] of the normalized inelastic cross section a nm = cr nm /a nm , 
where a nm denotes the average inelastic cross section in the energy interval AS (that is 
small compared to classical energy scales). This expectation is clearly confirmed by our 
numerical data presented in Fig. ^p)- 

The central figure of merit to identify Ericson fluctuations is the energy 
autocorrelation function 

C nm (e) = / dS(a nm (£ + e) - a nm )(a nm (£) - a nm ) . (8) 
J as 

A least-square fit of the numerically obtained autocorrelation as depicted in Fig. [5] shows 
that it perfectly matches a Lorentzian 

r 2 

C nm (e) oc £2 + v2 , (9) 

with mean resonance width T = 3.7-10~ 3 j[§J which is one order of magnitude larger than 
the mean level spacing A w 5 ■ 10~ 4 J, directly extracted from our numerical data. This 
is in perfect agreement with Ericson's scenario of overlapping resonances, and can be 
underpinned by a semiclassical picture [29] : The autocorrelation ([9]) can be interpreted 

§ We verified that a preparation in the regular regime (i.e., u -C 1 or u 3> 1), C nm (e) shows significant 
deviations from a Lorentzian. 
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Figure 5. (color online). Autocorrelation function C nm (e) (black o) calculated from 
the inelastic scattering signal shown in Fig. ,4b). The curve nicely matches a Lorentzian 
fit (blue). Inset: Semiclassical decay constant f3 versus the mean resonance width T 
(black □). The data points correspond to different values of the coupling constant 
7. The semiclassical result is obtained after averaging over several initial conditions 
in the chaotic regime. The predicted correspondence (/3 = T) is confirmed by the fit 
P = 0.92 T (red). 



as the squared Fourier transform of the survival probability P(t) of the probe particle 
to stay a given time t in the scattering region, i.e. on the TB site j = 0, hence with 
P(t) = |co(t)| 2 . That latter quantity is evaluated by direct solution of the classical 
evolution equations derived from (|3j [5]) (with initial conditions -P(O) = 1 and the GP 
system prepared at an energy corresponding to E m ), and exhibits an exponential decay 
P(t) oc e _/3i . (3 thus determines the width of the (classical) autocorrelation function 
C c i as (e), by virtue of 

2 



a 



clas(£) 



dt P(t) exp(iet) 



(10) 



which implies an average over all outgoing probe energies e n . The inset of Fig. [5] 
demonstrates a perfect match of this semiclassically extracted quantity f3 with the width 
T of the autocorrelation extracted from the quantum mechanical cross section, and thus 
provides an independent, semiclassical proof of the Ericson scenario in the present many- 
particle problem. 
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6. Conclusions 

In the light of recent BEC experiments, the proposed scattering setup represents an 
experimentally feasible and robust way to non-destructively probe the condensate. For 
example, a typical lattice depth of V = 10 E R recoil energies yields a tunneling strength 
K = 0.022 E R (cf. Eq.(20) of [3J). The exemplary case studied above (u = 5, N = 38) 
leads to a mean level spacing of A = 2.5 • 10 -3 E-&. This, in turn, corresponds to an 
energy resolution of A/E ~ 10 -3 needed to resolve single scattering channels, i.e. the 
final state of the system. As chaotic spectral statistics extend over a large parameter 
regime [IHl ED] one could also choose larger lattices (L > 3) and/or fewer bosons to 
increase A, and still obtain comparable results. Let us also stress that, in the proposed 
measurement, it is not necessary to determine the boson number on the lattice: As long 
as the detector can resolve the energy of the probe particle, the target properties can be 
inferred from the scattering signal, without knowledge of the exact number of bosons in 
the system. 

As expected, the observed chaotic scattering signals exhibit an inherent sensitivity 
on, e.g., the probe particle's energy. To obtain a reliable characterization of the target, 
we employed a statistical analysis based mainly on the autocorrelation function C(e) 
and the probability distribution P(ay) of the inelastic cross section. On the other 
hand, the robustness of the obtained results against unavoidable fluctuations in the 
experimental control parameters is corroborated by averaging our results over Am = 30 
channels (see, e.g., Fig. |4]). This averaging, together with the statistical independence 
of the S'-matrix elements, implies that our observations are robust against uncertainties, 
e.g., in the initially prepared state \E m ) of the condensate or in the particle number N, 
as well as in the resolution of the scattering channels. For a typical laser wavelength 
of 1064 ran, N = 38 Cs atoms, and u = 5, this average corresponds to a residual 
temperature of ~ 5 nK what is readily achievable with state-of-the-art experiments 
that reach temperatures as low as 0.35 nK [31]. 

Hence, measurements of the partial inelastic cross section can identify an 
unambiguous case of Ericson fluctuations, what yields robust information on the many- 
body spectrum. In contrast to compound nuclear reactions, here the latter is under 
perfect control, through the accurate experimental manipulation of the underlying Bose- 
Hubbard Hamiltonian, via the control parameter u. This would allow to experimentally 
investigate the fate of the Ericson fluctuations at the transition from the chaotic to the 
regular regime. By virtue of this control, our proposed setup could also add to the 
recent debate on complex many-body scattering (321 E31 HH] • Finally, the possibility to 
record single scattering channels could help the understanding of a related phenomenon: 
Summing up an increasing number of contributing channels would then allow to study 
in a controlled way the crossover from the Ericson regime to the multi-channel regime 
of universal conductance fluctuations. 
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